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Abstract 

Starting with the relative entropy based on a previously proposed en- 
tropy function S q [p] = J dx p(x)(—\np(x)) q , we find the corresponding 
Fisher's information measure. After function redefinition we then max- 
imize the Fisher information measure with respect to the new function 
and obtain a differential operator that reduces to a space coordinate sec- 
ond derivative in the q — * 1 limit. We then propose a simple differential 
equation for anomalous diffusion and show that its solutions are a gener- 
alization of the functions in the Barenblatt-Pattle solution. We find that 
the mean squared displacement, up to a (/-dependent constant, has a time 
dependence according to < x 2 >~ K 1 ^ q t 1 ^ q , where the parameter q takes 
values q — (superdiffusion) and q = ^Z-i (subdiffusion), Vn > f . 
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1 Introduction 

During the past two decades the Fisher information measure 1] has become, due 
mainly to the work of Frieden and coworkers, to play a fundamental role from 
where the well known covariant field equations of physics can be extracted [2]-[4] . 
Although the Fisher information measure is obtained from the relative entropy, 
it plays a different role than the Shannon's information measure. As physics 
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is concerned, in statistical mechanics maximization of Shannon's entropy under 
some constraints gives us the probability distribution function and therefore 
the partition function for a particular system, while in the Fisher information 
measure after defining a new function 4>(p) in terms of the probability amplitudes 
p(x), one obtains a 'kinetic' term in terms of (j> leading to the interpretation of 
the Fisher information as a kind of lagrangian. In particular, in Refs. 
it has been shown that the familiar relativistic wave equations are obtained 
maximizing the expressions I — J, where I is the Fisher measure in term of 
the new functions and J is obtained from the Fourier transform of I and use 
of the energy constraint. Curiously, this approach works very well for the case 
of covariant formulations. For the nonconvariant case, as it is the case of the 
Schrodinger and the diffusion equations, the Fisher measure defines the space 
coordinates differential operator and thus the time part has to be introduced 
by using dimensional analysis. 

This paper is organized as follows. In Section ((2|) we review the relation of 
the Fisher measure and the Kullback entropy, and briefly explain how the defi- 
nition of the new function as a function of p(x) leads to the interpretation of I 
as a lagrangian. In Section ([3]) we introduce the relative entropy for our entropy 
function and we obtain, following the method in Ref [5]), the corresponding 
Fisher measure, and then after function redefinition we find coordinate differen- 
tial operator. In this section we also obtain the eigenfunctions of this differential 
operator. In Section we show that this formalism the corresponding leads to 
an alternative differential equation for anomalous diffusion. In Section ([5]) we 
discuss our results. 



2 Fisher information measure as a lagrangian 

The Fisher information measure is defined through the integral 

/ dp(x 



I = dx 



\ dx 



(1) 



where p(x) denotes the probability distribution. In Ref. [5] it was shown that 
/ can be obtained by considering a small global shift A of p(x) in the Kullback 
information as follows 



(2) 



1(A) = J dxp(x) ln\p(x)/p(x + A)], 
and expanding in series in powers of the shift A, leading to the equation 



7(A) 



dx 



dp 

dx 



(gr Ma* 

V dx 2 ) 2 



(3) 



where one can recognize the Fisher measure in the quadratic term. In 1943, 
Fisher introduced [6] in Equation ([!]) a new function <^(x), called probability 
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2 



amplitudes, in terms of p(x) through the relation p(x) — 4> 2 {x) leading to the 
simpler equation^] 

1 = 4 fdx^ d t (X \ (4) 
J ax ax 

which is the starting point of Frieden's approach [3] to derive the equations that 
describe the behavior of physical systems. Introducing space-time coordinates, 
allowing two real components cj>t and </>2, and after defining $ = <fri + i<p2 
Equation ([4]) was promoted to a covariant expression \A\ 

I = 4 J d A xd v $*d v $, (5) 

such that performing the Fourier transform of I give us the energy constraint and 
thus the mass term m 2 $*$, leading to the well known Klein-Gordon equation. 



3 Fisher information based on fractional calcu- 
lus inspired entropy 

In a previous article [7] , based on the definition of the Riemann-Liouville frac- 
tional derivative, we proposed a new entropy function defined as 

5 g [p]= t lim i |^ 00 ^r 1 E e " lnPs )' ( 6 ) 
where < q < 1, and the fractional derivative is given by 

aA*= (^) n (-£>?""/(*)), (7) 
with n G N and n > q. The equation 

defines the fractional integral operator. Equation © leads to the entropy func- 
tion 

i 

which is a concave function and satisfies the Lesche and thermodynamic stability 
criteria. As explained in [7], the motivation to propose Equation (|6|) as a defi- 
nition of a new entropy function is based on the observation that the Shannon 
entropy and Tsallis entropy can be defined according to [8] 

S = lim — pj l 

i 

S q = t lim£><X>r*, (10) 

i 

2 Fisher used the symbol q(x) instead of <j>(x), which we will not use to avoid confusion with 
the parameter q. 
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where the operator 



i ri^df dt 

D\ = t-^—« . 

q 1-9 



(11) 



is called the Jackson q-derivative \9 . With use of Equation © for two prob- 
ability distributions p = (pi,...) and P = (Pi...) we define the corresponding 
relative entropy as 

ff(p||P)=5>(-ln^y. (12) 

The proof that H(p\\q) is positive follows along the same lines as the q = 1 case. 
Provided that '^2 i p% = 1 and f(x) is a convex function, the Jensen inequality 



(13) 



and the fact that f(x) = (— \nx) q is a convex function allow us to write for 
Equation (fl2|) . 



H(p\\P) > 



(- ta ?^)' 

(-m£P) 9 



(14) 



= 0. 



Now, to extract the Fisher information measure for q ^ 1, we consider a distri- 
bution function p(x) and a small shift A such that Equation (TT2"]) becomes 



P(p(ir)||p(.T + A)) = y dxp(x) |^-ln 
Expanding ln(p(x + A) up to second order in A, 



p {x + A) 
p(x) 



(15) 



P(p(x)||p(a; + A)) = - / dxp(x) 



P P A 



(r 


(t_ 


p l2 \ 


i 


n 




A P 


P 2 J 







where p' = -r-. In order to obtain the Fisher information it is enough to keep 
the two lowest order terms, 



H(p(x)\\p(x + A)) = — J dx ((«) J^A? - (?) 



p"p'1~ 1 \ A q+ 



r.q-1 



such that, up to a constant, we define 



/„ = dx 



p'i +1 {x) 
pi(x) 



(16) 
(17) 
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In our case, after introducing the new function 4>{x) 
becomes 



I q = (q + / dx 



which reduces to Equation ([4]) when q 
leads to the equation 

(d4>(x)\ q ~ l d 2 cj)(x 



dx J \dx 
■ 1. Setting 8I q 



1 (x), Equation (JT7J) 
(18) 

in Equation ([181) . 



\ dx 



dx 2 



= 0. 



(19) 



Although solutions to Equation (119)) are not different than the q = 1 case, 
solutions to the corresponding eigenvalue problem are different. Let us consider 
the equation 



(21) 

Taking the elementary 
, where solutions 



V dx J dx 2 - M{X) - 
Solutions to Equation (|2"D)) are given by the function 

cj>{x) = A [1 + (1 - 9)70;]" , 
where A, 7 and w are constants to be determined, 
derivatives we find that uj = — f~ an d 7 = ( 2(9+ 1)^- 1 

g + l 

with 7 > restrict the values of q such that (—1) 2q = 1. It is interesting to 
remark that functions in Equation (f2"Tj) are of the type of solutions found for the 
probability distribution in the non-extensive formulation of statistical mechanics 
[1011. The constant .4 is found by normalizing the function in Equation ( 21) 
using the definition of the T function, in a very similar way as done to calculate 
the partition function of non-extensive statistical mechanics |12j [13j , 



r(z) 



"dt x > , z > 0, 



(22) 



leading to the value A = 7(1 + q). 

Taking the limit q — -> 1 the function in Equation (|21|) becomes the expected 
exponential solution 4>{x;q = 1) = A 1 / 2 exp(— A 1 / 2 ^). 



4 Anomalous diffusion 

In the last twenty years, anomalous diffusion has been observed in a variety of 
interesting physical systems, for example breakable micelles dissolved in salted 
water [14], measurements of cardiac beat to beat intervals in healthy individu- 
als |15j . chaotic transport in a two-dimensional fluid flow in a rotating annulus 
[16] , laser cooling of atoms to very low temperatures |17j . and more recently 
several interesting experiments with living cells in which a transition from sub- 
diffusive to superdiffusive behavior is observed [IS] 02]. On the other hand, 
several mathematical models have been proposed to study different aspects of 
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anomalous diffusion. Some of these models are based on a linear differential 
equation for diffusion on fractals [20] . linear and non-linear fractional differen- 
tial equations [H]-[27], a study of particle chaotic dynamics along a stochastic 
web |28j . a study of super-diffusion in a hamiltonian system in the context of 
a continuous-time random walk approach to Levy flights [53], nonlinear differ- 
ential equations [3^-[3S] and some calculations based on Tsallis non-extensive 
statistical mechanics [36]-[40]. 

In our case, the simplest model we can propose is based on the differential 
equation 

where K is the diffusion constant. The solution to the eigenvalue problem in 
the previous section guides us to propose the function 

W(x,t) = ^(x,tr, (24) 



where 

. . I" (1 - q) 



1 



C 



(25) 



and A, A, C, (3, 7 and to are constants to be determined, which in the q — > I 
limit must satisfy: A — > 1/2, j3 — ► 1 and 7 — > 2. By performing the elementary 
derivatives, comparing powers in ^(x, t), x and i and the constants in both sides 
of Equation (f23f we find that 

" " -T^ < 26 > 

2q 2 

A = i (29) 



/3 = & (28) 



2g 

x _ q _ 2^(g + i; 



1 



(30) 



Similarly to the previous section there is a restriction on the values of q, due 
to the fact that in order to have the correct q — * 1 limit the constants A and 
C have to be positive leading to the set of allowed values of q = 2n+i w -h ere 
n = 1,2,3.... Since the solution W(x, t) must have dimension [W] = these 

constants have dimensions: [K] = £JI ^-, [A] = and [C] = sn §-. In Fig. 
1 we show the function W(x,t) for the case of q = 5/7 and t = 0.05,0.09 
and 0.2, and in Fig. 2 the function W(x, t = 0.5) for the superdiffusive cases 
q = 5/7, q = 19/21 and the normal diffusive case q = l. 
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As a next step we find an additional relation between the constants A and 
C by normalizing the solution W(x, t) 



1 



(1-9) 



C ifl 



= i. 



(31) 



Since q takes values such that (— l) 7 = 1 we can change the limits from (— oo, oo) 
to (0, oo), such that with use of the integral representation of the T function we 
obtain 



.4 



1 



C 

- ' (-0)) 



1/7 



.9(<7), 



(32) 



where the function g(q) = 2 r(i/7)r(-J-i/ 7 ) ■ Dennin g Q = and b Y usin S 

m 

lim T<yQ+ Il e-^^Q - 1 (33) 
T(Q) 6 

with Equation (l30|) we find that the constant A has the correct normal diffusion 
limit A — > \J~j^j< as g — > 1. The constant C is given by 

C* = (2^)^(1 + g) 7/2 (l - g) 1/2w .g 7/2w (q). 
The mean squared displacement < .x 2 > is obtained by solving 



(34) 



< x >= 



.4 

7^ 



dx 



1 



(1 - g) x 7 



C f 3 



leading, up to a g-dependent constant, to 

< x 2 >~ K 1 / q t 1 /\ 0<g<l. 



(35) 



(36) 



Since < q < 1, our model seems to describe exclusively superdiffusive pro- 
cesses. This apparent drawback is due to the fact that we took in Equation 
^ the particular value n — 1. A general expression, Vn > for the entropy 
function S q [p) reads 



-t In pi 



which leads to the same mathematical function 



(37) 



(38) 



but now the parameter q takes values in the set n — 1 < q < n. The solution 
for q > 1 reads 



W'(x,t) 



A[_ 

= 



1 - 



(9 - 1) x< 



C> ifi 



for(l-^^)<0, g>l, 



(39) 
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with the constants A' and C given by 



.4' 

c 



7A 



u,iK(j - 1) 



A 



<?'(<?) 3 

X 

" (q-l) W V + V(<7)^, 



(40) 



with g'(q) = 2r(iJj)r(i+L) • ^ n or dcr to normalize the new solution we used the 
integral representation of the T-function in the complex plane 



T(z) 2tt 



iz * (-*)' 



dt. 



c 



(41) 



In this case, q > 1, the constant cj is positive and therefore the solution W'(x, t) 
becomes the usual normal diffusion solution in the q — ► I limit. The allowed 
values of q are q = 2—1 > an< ^ thus Equation (|36[) for an arbitrary positive integer 
n, includes also the case of subdiffusive behavior. Fig. 3 displays the function 
W'(x,t = 0.5) in Equation (J39J for the subdiffusive cases q = 7/5, q = 21/19 
and the normal diffusive case q = 1. 

It is interesting to remark that for the particular value 7 = 2 and arbi- 
trary values for the constants w, A and (3 the functions W(x, t) and W'(x, t) in 
Equations ([24]) and (|39[) are of the same type of the solutions of the nonlinear 
diffusion equation studied in [30] [32] [42] , called the Barenblatt-Pattle solution. 



5 Conclusions 

In this paper, based on a previously proposed entropy function S q [p] we first 
defined the corresponding relative entropy H(p\\P), and by considering the two 
probability distributions to differ by a small global shift A, we expanded the 
relative entropy with respect to this shift and extracted the Fisher information 
measure for q =/= 1. Then, following Frieden's idea we defined a new function 
4>(x) in terms of the probability p(x) that allowed us to reinterpret the Fisher 
information as a free action I q . Setting the variation of I q to 5I q = we find a 
differential operator that reduces to a second order derivative at q = 1. We found 
that the eigenfunctions of this differential operator are of the same type than 
the probability distribution found in the Tsallis formulation of noncxtcnsivc 
statistical mechanics. Using this differential operator we proposed a simple 
differential equation that coincides with the normal diffusion equation in the 
q — ► 1 limit. In our case, the solutions are such the exponent 7 in the x 
variable is a fraction 7 = taking values in the set 7 = „ 4 " . for q = 2 I L ~\ 
(superdiffusion), and 7 = 2 ^" 1 for q = 2Z-1 (subdiffusion),Vn = 1,2,3.... The 
work done in this paper is a first step to study anomalous diffusion with a rather 
simple and different mathematical model. In order to study anomalous diffusion 
in a more realistic way will require to introduce a time dependent source term 
(absorption) and external forces (drift) in a three dimensional version of this 
model, which we will attempt to address in future communications. 
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Fig. 1: The solution W(x,t) to Equation (j2"5)) for the case of q = 5/7 and 
t = 0.05,0.09 and 0.2. 



Fig. 2: The function W(x,t = 0.5) for the superdiffusive cases q = 5/7 
(thin solid \ine),q = 19/21 (dotted) and the normal diffusive case q = 1. 



Fig. 3: The function W'(x,t = 0.5) in Equation flU for the subdif- 
fusive cases q = 7/5 (thin solid line), q — 21/19 (dotted line) and the 
normal diffusive case 9=1. 
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